Simulating redox coupled proton transfer in cytochrome c oxidase: Looking for the proton bottleneck

AbstractGaining a detailed understanding of the molecular nature of the redox coupled proton transfer in cytochrome c oxidase (COX) is one of the challenges of modern biophysics. The present work addresses this by integrating approaches for simulations of proton transport (PTR) and electron transfer (ET). The resulting method converts the electrostatic energies of different charge configurations and reorganization energies to free-energy profiles for different PTR and ET pathways. This approach provides for the first time a tool to study the actual activation barriers and kinetics of different feasible PTR processes in the cycle of COX. Using this tool, we explore the PTR through the bottleneck water molecules. It is found that a stepwise PTR along this commonly assumed path leads to far too high barriers and is, thus, inconsistent with the observed kinetics. Furthermore, the simulated free-energy profile does not provide a simple gating mechanism. Fortunately, we obtain reasonable kinetics when we consider a PTR that involves a concerted transfer of protons to and from E286. Finally, semi-qualitative considerations of the forward and backward barriers point toward open questions about the actual gating process and offer a feasible pumping mechanism. Although further studies are clearly needed, we believe that our approach offers a general and effective tool for correlating the structure of COX with its function

Progress in the Prediction of pKa Values in Proteins

The pKa-cooperative aims to provide a forum for experimental and theoretical researchers interested in protein pKa values and protein electrostatics in general. The first round of the pKa-cooperative, which challenged computational labs to carry out blind predictions against pKas experimentally determined in the laboratory of Bertrand Garcia-Moreno, was completed and results discussed at the Telluride meeting (July 6–10, 2009). This article serves as an introduction to the reports submitted by the blind prediction participants that will be published in a special issue of PROTEINS: Structure, Function and Bioinformatics. Here, we briefly outline existing approaches for pKa calculations, emphasizing methods that were used by the participants in calculating the blind pKa values in the first round of the cooperative. We then point out some of the difficulties encountered by the participating groups in making their blind predictions, and finally try to provide some insights for future developments aimed at improving the accuracy of pKa calculations

Transition state theory can be used in studies of enzyme catalysis: lessons from simulations of tunnelling and dynamical effects in lipoxygenase and other systems

The idea that enzyme catalysis involves special factors such as coherent fluctuations, quantum mechanical tunnelling and non-equilibrium solvation (NES) effects has gained popularity in recent years. It has also been suggested that transition state theory (TST) cannot be used in studies of enzyme catalysis. The present work uses reliable state of the art simulation approaches to examine the above ideas. We start by demonstrating that we are able to simulate any of the present catalytic proposals using the empirical valence bond (EVB) potential energy surfaces, the dispersed polaron model and the quantized classical path (QCP) approach, as well as the approximate vibronic method. These approaches do not treat the catalytic effects by phenomenological treatments and thus can be considered as first principles approaches (at least their ability to compare enzymatic reaction to the corresponding solution reactions). This work will consider the lipoxygenase reaction, and to lesser extent other enzymes, for specific demonstration. It will be pointed out that our study of the lipoxygenase reaction reproduces the very large observed isotope effect and the observed rate constant while obtaining no catalytic contribution from nuclear quantum mechanical (NQM) effects. Furthermore, it will be clarified that our studies established that the NQM effect decreases rather than increases when the donor–acceptor distance is compressed. The consequences of these findings in terms of the temperature dependence of the kinetic isotope effect and in terms of different catalytic proposals will be discussed. This paper will also consider briefly the dynamical effects and conclude that such effects do not contribute in a significant way to enzyme catalysis. Furthermore, it will be pointed out that, in contrast to recent suggestions, NES effects are not dynamical effects and should therefore be part of the activation free energy rather than the transmission factor. In view of findings of the present work and our earlier works, it seems that TST provides a quantitative tool for studies of enzyme catalysis and that the key open questions are related to the nature of the factors that lead to transition state stabilization